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We calculate a nonequilibrium steady state of a quantum XX chain in the presence of dephasing 
and driving due to baths at chain ends. The obtained state is exact in the limit of weak driving while 
the expressions for one- and two-point correlations are exact for an arbitrary driving strength. In 
the steady state the magnetization profile and the spin current display diffusive behavior. Spin-spin 
correlation function on the other hand has long-range correlations which though decay to zero in 
either the thermodynamical limit or for equilibrium driving. At zero dephasing a nonequilibrium 
phase transition occurs from a ballistic transport having short-range correlations to a diffusive 
transport with long-range correlations. 

PACS numbers: 75.10.Pq, 05.60.Gg, 05.70.Ln, 03.65.Yz 
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Introduction- The theory of equilibrium processes is 
well developed with known equilibrium distribution func- 
tions. Nonequilibrium physics is on the other hand much 
more difficult to treat. One of the simpler nonequilib- 
rium situations is the steady state setting to which the 
system relaxes after a long time under nonequilibrium 
driving. Still, it is not even known if a generic distribu- 
tion function for a nonequilibrium steady state (NESS) 
exists. Particulary important aspect of NESS is the ques- 
tion of transport, for instance, it is not fully understood 
under what conditions does one have a diffusive trans- 
port [l|. There is also a large disparity between classi- 
cal and quantum systems. For many classical models, 
like various exclusion processes [4I, an explicit solution 
is known, greatly facilitating our understanding of inter- 
esting phenomena that can occur in NESS. Physics of 
quantum NESS is on the other hand still a rather un- 
charted territory with only a few exact solutions. One 
such class of exactly solvable systems are those described 
by master equations quadratic in fermionic operators. 
NESS in doubly-infinite XY chain has been studied in [3|. 
Quadratic open systems can in fact be exactly diagonal- 
ized in operator space [J| making it possible to study 
nonequilibrium phase transitions [5|. XX chain interact- 
ing with baths has been analytically studied also in [6|. 
Time evolution of quadratic open systems can as well be 
calculated efficiently in terms of matrix product opera- 
tors [7| , similarly as for closed systems [8J . Tight-binding 
(ie., XX) model in the presence of the environment has 
been studied in j9|. Superoperator corresponding to the 
master equation has been diagonalized and found that 
there is a crossover from nondifFusive to diffusive behavior 
as the chain size is increased. Diffusion in an XX model 
with classical noise has been considered in [10|. In a re- 
cent work [111] quantum exclusion processes, described by 
the master equation with nearest-neighbor Lindblad op- 
erators for stochastic jumps as well as with coherent uni- 
tary XX part, have been numerically studied. For some 
parameters similar long-range correlations are obtained 
in the NESS as here. 

In the present work we are going to provide an ex- 



plicit solution for the NESS of a quantum model, show- 
ing that the system exhibits diffusive transport and has 
long-range correlations. Furthermore, as a parameter of 
the model is varied a nonequilibrium phase transition oc- 
curs. To our knowledge this is the first quantum model 
for which one is able to analytically show diffusive be- 
havior. 

The system studied is a one-dimensional XX chain of 
spin-(l/2) particles, 



^ = EK^^i + ^J-J+i)- 
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To induce a nonequilibrium situation we couple the sys- 
tem at both ends to baths at different potentials. In 
addition, each spin is also exposed to a depha sing . The 
same model has been numerically studied in 13|. The 
evolution of system's density matrix is governed by the 
master equation of the Lindblad form IJ] , 
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Nonunitary terms C^'^^'^ and ^'^'^p'^ are expressed in terms 
of Lindblad operators as Ylkyl-^kPiLki + lLk, pL^^, 
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a — a'^ ±ia'^. Two parameters for baths are the cou- 
pling strength F and the driving strength fi. Dephas- 
ing acts independently on each spin separately, C^^p^ — 

^(T^. It is diagonal in the Pauli ba- 



^ j^ £ , with £ ■ having only one Lindblad op- 

, T dcph 

erator L^ '■ — 

sis, £fP''(af^) = -27ar^ while £r^"(a?) = and 



/>dcph 



Ij) = 0. Dephasing with strength 7 therefore 



causes an exponential decay of the off-diagonal elements 
in the diagonal basis of a'' . 

To find the NESS of the master equation ^ we are 
going to write p as a sum over 4" different products of 
Pauli matrices and solve for unknown coefficients by de- 
manding stationarity C(^p) = 0. Using Jordan- Wigner 
transformation we can express the whole system in terms 
of spinless fermions. While the H ([T]) and /^t"^*'' ([3]) are 
quadratic in fermionic operators, £deph jg ^^^^ Because 
it involves a product of two a"^ it is quartic and there- 
fore the recently introduced [J, |5| analytical solution for 
open quadratic systems can not be used. Note that all 
quadratic systems are ballistic. 

Let us start with an unnormalized ansatz, 
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with the individual terms being 
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where j^ = 2(tT^cr^_^-^ — cr^cr^+i) is the operator of local 
spin current, while the other two terms are 
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We shall show that using an appropriate coefficients in 
the above ansatz, acting with C on such p gives zero up-to 
0{p^), ie., that the above ansatz is the correct solution 
up-to 3rd order in the driving strength p. Powers of p 
written in front of various coefficients in the ansatz sug- 
gest a scaling of the corresponding coefficient in the solu- 
tion, that is, post festum we shall see that indeed pA ex p, 
p'^C oc p"^, etc.. Note that for equilibrium driving, ^ = 0, 
the solution is p '-^ 1, which can be interpreted as an in- 
finite temperature state, although one must be careful 
speaking about the temperature because integrable sys- 
tems in general do not thermalize, even when coupled to 
reservoirs J12l] . 

Demanding that the expansion coefficients in front of 
all operators in C{p) vanish one obtains a set of equations. 
Coefficients in front of two boundary operators af „ give 

- 6 - T/i - Tai = 0, b + Tp - Tan = 0. (9) 

These two equations are exact to all orders in p because 
one can get a^ „ term in C{p) only from terms of the form 



1, al „ or ji^n in P, which are all included in our ansatz 
Q! That is', tT[alnC{p)] y^ only if p € {l,<Tf^„, ji,„}. 
As we shall see, all equations we are going to use are, in 
fact, exact to all orders in p\ Equations ([9]) give relation 
a„ — — oi and oi = —p— y- From coefficients in front of 
a^j in the bulk (j 7^ f , n) one gets that the expectation of 
the local current jk is independent of the site fc, ie. bk — 
b, which has already been explicitly taken into account 
in our ansatz. 

Equations which we get by demanding that the coeffi- 
cients in front of ji are zero are as follows: 
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with T^^'^^ = 27 -f TSi^i + TSi+k-i,n, while lif' is the 
efficient in the expansion of p in front of the operator 



^3 - '^jf^J + l'^J+2 
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At the boundaries, 



terms from operators having site indices out of chain 
range 1 < j < n are absent. These equations pUj) are 
exact because we can get jk term only from jk , cr^ or 

f^fe-i^fe^fe+i + ^l-i^l^l+i- Before solving Eq.dTU]) we 

(3) 
are going to show that ft,) ' are actually zero. 

To this end let us in addition to the expansion co- 
efficients which are explicitly included in the ansatz. 



also denote by h, the expansion coefficient in front 



of the operator H- ' — a^a" 
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(TjCr^j+i ■ ■ ■ o-|_|_j._2 0'J+fe_i- B^' is a generalization of the 

(k) 

spin current while if is the hopping operator to the 
{k — l)th neighbor. Note that the rotational symmetry 

of C around z-axis imposes the value of the two signs 

(2) 
in both operators. For b\. we already know that it is 

site-independent and equal to 6, in addition we define 

(1) (k) 

b). = 0. We are going to show that all h^ ' as well as 
6 J- (apart from h^ ) are zero. Demanding that the co- 
efficients in front of i7- in £(p) are zero, we get for 
k = 2, . . . ,n the equations 
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Equations from the coefficients in front of B- are, on 
the other hand, for fc = 3, . . . , n 

{h^+'^-h^^) + {h!f-'^-hf;^^)+i^hf = 0. (12) 

With exact equations pT|) and p2p we have a closed set 
of exactly as many equations as there are unknown h- 

(k) 

and b- . Homogeneous set is nonsingular with the only 

(k) (k) 

solution being trivial h, = and b, =0. This means 

(k) 

that NESS p does not contain any B^- ' terms with fc > 2 

(k) 

nor any H^j ' terms with fc > 2. 



Going now back to the equations obtained from cur- 
rents ([To]), they can be solved exactly. Together with two 
equations © they give the solution for b and all Gj, 
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O-n-l 



-6(-+r + 2(n-2)7)-Ai 



a„ = -6(- + 2r + 2(n-l)7)-M. (14) 

This solution is exact to all orders in fi for any n, ir- 
respective of other terms present in the expansion of p 
not included in our ansatz. Before proceeding to solve 
for C, D and F let us for a moment discuss its physi- 
cal significance. First, observe that the current and the 
magnetization expectation values in the NESS are sim- 
ply equal to the corresponding coefficients, (a^) — aj, 
(j) = 2b. Magnetization has a linear profile and the cur- 
rent scales as j ~ p,/n for large n as long as 7 ^ 0. 
Such behavior is typical for a system with normal (dif- 
fusive) transport. Spin conductivity is k = I/7 and as 
one decreases the dephasing strength 7 the current di- 
verges as ~ 1/7 [l3|. Also interesting, the limit of weak 
coupling to the bath, F — > 0, is at fixed n and zero de- 
phasing singular. That is, the current is proportional to 
F as J ^ fjT while the magnetization at the boundaries 
is proportional to /iF^. In particular, without dephasing 
for 7 = 0, transport is ballistic with b = — F/i/(l + F^), 
ai = -F2^/(1+F2), a„ = F2^/(1+F2), while aj^i,„ = 0, 
the same as in the recent result in 6]. 

To avoid cumbersome expressions for the 2nd order 
terms we from now on set 7 = F = 1. This has no qual- 
itatively relevant consequences. In this case the solution 
(I14I13P can be simply written as. 



fJ- 



1 - k. 



71+1 ■' n + 1 

Expansion coefficients in front of <j\cr\ give (using (|15l) ). 



-p., kj = 2j-Sj^i+Sj,n. 
(15) 



36^ + di + Ci.2 =0, for fc = 2 

dn-2 ~ C'l.s ~ 0, for fc = 3 

d„_3 - Ci,4 = 0, for fc = 4 

di — C'l n — 0, for k — n. 



Interior nearest-neighbor terms a'^a'^^j^^ give (non nearest- 
neighbor are zero) 

26^ + dj -f dn-j = 0, for j = 1, . . . , n/2. (17) 



Equations (|16ll7p are again exact because cr^cr^ term can 
only be obtained from terms of the same form or from 
jk<J^j or (7|. Coefficients in front of (j\ik, on the other 
hand, give 

f -2b^ + 4di - Ci,2 + Ci,3 = 0, for fc = 2, 
2f - 2b^ + Ml - Ci,k + Ci^k+i =0, 2<k<n-l, 
2/ - 26^ -H 5di - Ci,„_i + Ci^„ = 0, k^n-1. (18) 

Equations are again exact for the same reason as be- 
fore. The set of equations (|16ll7ll8p can be solved, giv- 
ing explicit expressions for di, f and Ci^k- Using these 
in equations obtained for (T2Jk enables one to obtained 
C2.J, then from a'^jk we obtain Csj and so on. At the 
end the solution is. 
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C,j = (/c,/c„+i_j + (n + 1) • 5j- ,+1), (19) 

n 

where kj is the same is in aj (|15p . The solution is exact 
to all orders in p,. One can check that the ansatz with 
coefficients (fTO)) and ((T5|) results in the C{p) having only 
terms of order 0{^^). One could systematically calculate 
also 3rd and higher order terms, however, because main 
physical picture already emerges from the 2nd order term 
we are going to stop here (due to algebraic structure all 
nonzero operators in C{p) are in fact products of three 
terms, each being either cr^ or jk] more, 3rd order term 
in the ansatz would itself be a sum of such products of 
three operators, each being either jk or a'-). 

Because coefficients in the ansatz are exact to all or- 
ders in p, the corresponding single- or two-point expecta- 
tion values are also exact. Two-point connected correla- 
tion function of the magnetization, C{i,j) = (a^a^) ~ 
(cr^)(crp, is symmetric across the diagonal and is for 
j > i negative and equal to Cij p^ . On the diag- 
onal we have C{i,i) = 1 — af. If we are away from 
the diagonal, j j^ i + 1 and i,j 7^ l,ri,, the correlation 
function is, in terms of scaled variables, simply equal to 
Cix = i/{n + 1), 2/ = j/{n + 1)) = -^x{l ~ y). 

To independently verify our exact analytical solution 
for magnetization profile and correlation function we 
have compared it to the numerical simulation. Using the 
tDMRG 17[ we can calculate numerically exact expecta- 
tion values in the NESS for chains with n ~ 100 spins. In 
Fig. [I] one can indeed see that the numerical and analyt- 
ical solution agree perfectly also for very strong driving 
/i = 0.9. 
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exists a deeper connection with classical statistical mod- 
els remains to be seen |ll| . All these results suggest that 
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FIG. 1: Comparison between the numerical calculation (sym- 
bols) for n = 32 and strong driving ^ = 0.9 and the analytical 
solution (Ea ll5|[T9l full line). Top frame shows the magneti- 
zation profile and the bottom the correlation function. 



Correlation function indicates that there is long-range 
order present, which though goes to zero either in the 
thermodynamic limit n — > oo, or in the equilibrium 
limit /i — > 0. Long-range order is therefore a hnite size 
nonequilibrium phenomenon. Such correlations in the 
same XX model with dephasing have been numerically 
observed in [l3 |. There, a qualitatively same behavior 
has been obtained also for an anisotropic XXZ model 
with dephasing, perhaps suggesting that the form of the 
solution found here is more general and would approx- 
imately apply also for other systems. In particular, we 
have indications 15] that a similar behavior is obtained 
for pure XXZ model without dephasing, whose seemingly 
normal spin transport [16] in the gapped regime is sur- 
prising and needs further understanding. It appears that 
long-range correlations are a rather common feature of 
quantum NESS 0, [U, [H [H] . The form of the corre- 
lation function C{x,y) = — ^ '^' x{l — y) is the same as 



the form of the exact solution found here might have a 
more general validity. 

Let us for the end briefly discuss two-point correlations 
for the XX model without dephasing, 7 = 0. Calculation 
along the same lines as above results in the solution (for 
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together 



with the already obtained b = —^ and aj^i.„ = 0, while 
fli — — a„ — ~2' ^^ '^^^ ^^^ ^^^^ there are no long-range 
correlations present. Therefore, at 7 = a nonequilib- 
rium phase transition occurs from a ballistic transport 
with short-ranged correlations for 7 = to a diffusive 
regime with long-range correlations for 7^0. 

Conclusion- We have analytically calculated a 
nonequilibrium steady state of an open quantum spin 
chain to the 2nd order in the driving strength. One and 
two-point correlation functions are calculated exactly to 
all orders. The system, XX model with dephasing, shows 
diffusive transport with linear magnetization profile and 
current scaling as ~ 1/n. There are also long-range cor- 
relations present in the steady state. To our knowledge 
this is the first quantum system for which one is able to 
analytically show diffusive transport. At zero dephasing 
a nonequilibrium phase transition occurs from a short- 
ranged ballistic to a diffusive transport with long-range 
correlations. Furthermore, long-range correlations in the 
present model have the same form as in some classical 
statistical exclusion processes and there are some indica- 
tions that the same qualitative behavior is found also in 
other quantum systems. Therefore, while the exact so- 
lution might be possible only for this special model, the 
results might have more general (linear response) validity. 
It is hoped that the presented result will lead to better 
understanding of quantum transport in one-dimensional 
systems as well as of the conditions under which long- 
range correlations appear in nonequilibrium quantum 
states. This could bridge a gap between relatively well 
explored classical exclusion precesses, where many ana- 
lytical solutions are known, and quantum nonequilibrium 
steady states where almost no exact solutions are known. 
MZ is supported by the Program P 1-0044 and the Grant 
J 1-2208 of the Slovenian Research Agency. 



in many classical exclusion processes [2|. Whether there 
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